nb <- mgcv::gam(count ~ s(lon,lat, k = 10) +
s(lod, k = 5, bs = "cc")+
s(sst, k = 5)+
s(min_dist_owf, k = 10)+
s(min_dist_shipwreck, k = 10),
offset = log(active_tags),
data = detections_day,
family = "nb",
method = "REML")
library(glmmTMB)
zinb <- glmmTMB(
count ~ s(lon, lat, k = 10, bs = "tp") +
s(lod, k = 5, bs = "cc") +
s(sst, k = 5) +
s(min_dist_owf, k = 10) +
s(min_dist_shipwreck, k = 10) +
offset(log(active_tags)),
zi = ~ sst + min_dist_owf,
REML = TRUE,
family = nbinom2,
data = detections_day
)
zanb <- glmmTMB(
count ~ s(lon, lat, k = 10, bs = "tp") +
s(lod, k = 5, bs = "cc") +
s(sst, k = 5) +
s(min_dist_owf, k = 10) +
s(min_dist_shipwreck, k = 10) +
offset(log(active_tags)),
zi = ~ sst + min_dist_owf,
family = truncated_nbinom2,
REML = TRUE,
data = detections_day
)nb <- mgcv::gam(count ~ s(lon,lat, k = 10) +
s(lod, k = 5, bs = "cc")+
s(sst, k = 5)+
s(min_dist_owf, k = 10)+
s(min_dist_shipwreck, k = 10),
offset = log(active_tags),
data = detections_day,
family = "nb",
method = "REML")
library(glmmTMB)
zinb <- glmmTMB(
count ~ s(lon, lat, k = 10, bs = "tp") +
s(lod, k = 5, bs = "cc") +
s(sst, k = 5) +
s(min_dist_owf, k = 10) +
s(min_dist_shipwreck, k = 10) +
offset(log(active_tags)),
zi = ~ sst + min_dist_owf,
REML = TRUE,
family = nbinom2,
data = detections_day
)
zanb <- glmmTMB(
count ~ s(lon, lat, k = 10, bs = "tp") +
s(lod, k = 5, bs = "cc") +
s(sst, k = 5) +
s(min_dist_owf, k = 10) +
s(min_dist_shipwreck, k = 10) +
offset(log(active_tags)),
zi = ~ sst + min_dist_owf,
family = truncated_nbinom2,
REML = TRUE,
data = detections_day
)nb <- mgcv::gam(count ~ s(lon,lat, k = 10) +
s(lod, k = 5, bs = "cc")+
s(sst, k = 5)+
s(min_dist_owf, k = 10)+
s(min_dist_shipwreck, k = 10),
offset = log(active_tags),
data = detections_day,
family = "nb",
method = "REML")
library(glmmTMB)
zinb <- glmmTMB(
count ~ s(lon, lat, k = 10, bs = "tp") +
s(lod, k = 5, bs = "cc") +
s(sst, k = 5) +
s(min_dist_owf, k = 10) +
s(min_dist_shipwreck, k = 10) +
offset(log(active_tags)),
zi = ~ sst + min_dist_owf,
REML = TRUE,
family = nbinom2,
data = detections_day
)
zanb <- glmmTMB(
count ~ s(lon, lat, k = 10, bs = "tp") +
s(lod, k = 5, bs = "cc") +
s(sst, k = 5) +
s(min_dist_owf, k = 10) +
s(min_dist_shipwreck, k = 10) +
offset(log(active_tags)),
zi = ~ sst + min_dist_owf,
family = truncated_nbinom2,
REML = TRUE,
data = detections_day
)